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Abstract 



' We systematically investigate the embedding of U{l)x Froggatt-Nielsen models in 

(four-dimensional) local supersymmetry. We restrict ourselves to models with a 
single flavon field. We do not impose a discrete symmetry by hand, e.g. imparity, 
CN ' baryon-parity or lepton-parity. Thus we determine the order of magnitude of the 

in 
o . 

, scenario. We then scrutinize whether the predicted coupling constants are in accord 

with weak or GUT scale constraints. Many models turn out to be incompatible. 

p_; 1 Introduction 

43 ■ 

The problem of the fermionic mass spectrum is unresolved. Within the Standard Model 
(SM), the masses are due to the Yukawa interactions of the Higgs scalar. A theory of 
fermion masses is thus a theory of the origin of the Yukawa couplings, a problem not dealt 
with by the SM. This persists when extending the SM to supersymmetry. In fact the 
problem becomes worse. Including all terms in the superpotential which are allowed by 
the gauge group, there are 45 unknown extra .R-parity violating Yukawa couplings 
beyond the 27 of the SM. A theory addressing the fermionic mass spectrum should 
explain all parameters in the superpotential, i.e. the Jfi p coupling constants as well as the 
bilinear superpotential terms (the /i-problem) and possible higher- dimensional operators. 
There is extensive experimental data on the conservation of baryon- and lepton-number 
and thus bounds on the fi p coupling constants. It is the purpose of this paper to show 
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how this can be utilized to constrain the supersymmetric flavour theory, possibly pointing 
towards a solution of the fermionic mass problem. 

The problem of the fermionic mass spectrum has been extensively discussed in the 
literature. In Refs. [1-16] and Refs. [17-21], a certain class of (mostly) supersymmetric 
Froggatt-Nielsen models j22] is presented. The models in Refs. [1-16] all feature: 

1. A single local U(l)x symmetry. The gauge charges, X, of the superfields are 
generation-dependent . 

2. This U(l)x gets broken spontaneously by a single scalar flavon field. 

3. The mixed U(l)x anomalies are canceled by the Green-Schwarz mechanism |2*3*] . 

4. The fermionic mass spectrum at the grand unified (GUT)-scale satisfies |24| I25j . 

m e : : m T 

m T : nib 
m d :m s : m b 

m b : m t 
m u :m c : m t 
m t 

jjCKM 

Here tan j3 is the ratio of the vacuum expectation values (VEVs) of the two neutral 
Higgs scalars, tan/3 = (H u ) / (H v ); JJ CKM is the Cabibbo-Kobayashi-Maskawa ma- 
trix and e ~ 0.22 is the Wolfenstein parameter, i.e. the sine of the Cabibbo angle. 
Since the expansion parameter is rather large (e and e 2 are of the same order of 
magnitude) the powers of e are only approximate. 

5. The models are all in four-dimensional space-time. 

The models in Refs. [1-9] do not contain right-handed neutrinos, the models in Refs. [10- 
16] do. We have tried to present an exhaustive list of the existing models fulfilling the 
four points stated above. One model in Ref. [16] and the models in Refs. [17-21] fulfill 
the first three points but deviate in their fermionic mass spectrum. 1 

In this paper, we present a systematic embedding of Froggatt-Nielsen models in su- 
persymmetry, including a detailed discussion of the Kahlerpotential. We then investigate 
whether the U(l)x charge assignments of the aforesaid models give rise to an acceptable 
low-energy phenomenology. We hereby make the following assumptions: 

1 We do neither consider a gauged U(l) family-dependent i?-symmetry as in Ref. [2H] nor a discrete 
gauge symmetry as in Ref. |27j . 



~ e 4or5 :e 2 :l, (1.1) 

~ 1, (1.2) 

~ e 4 : e 2 : 1, (1.3) 

~ e ' 1 ' 2 or 3 cot (3, (1.4) 

~ e 8 :e 4 :l, (1.5) 

~ (H U ), (1-6) 
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1. The models are supersymmetric. 

2. Baryon- [2E] and lepton-parity (the latter one being anomalous) are not imposed 
by hand, so a priori all superpotential terms that are allowed by the symmetries 
are assumed to indeed exist; we shall hence include (for an introduction see e.g. 
Refs. [23 Hill inH E2|) and also higher dimensional operators. 2 

3. Supersymmetry is local and the breaking of supersymmetry is mediated by gravity 
[SSI IHZl BHI IHSj; hence the gravitino mass M 3 / 2 is of the same order as the masses 
of the squarks and sleptons, and the Kim-Nilles [10] or the Giudice-Masiero jH] 
mechanism solves the //-problem [this is consistent with the existence of U(l)x, see 
Eq. (EIZI)]- 

When the U(l)x is broken (c.f. Section|2j), the i?-parity violating minimal supersymmet- 
ric standard model is generated. In general, we shall allow for three right-chiral neutrino 
superfields and we denote the model by fi p -MSSM+N l . Its most general renormalizable 
superpotential is given by 



£ ab §X y G (U) Qi jjj 



+ e ab 5 xy G { % Q 
+ e 



H\ D\ 



ab L\ 



H v b E3 



+ e ab ii H u a H' 



v 



T U 



+ e ab L\ 



+ 



Tij N { Ni 



+ - e ab A ijk L\ L\ 



+ e ab 5*y A' ijk 



+ - e xyz A"ijk U\D\ D\ 



+ e 



ab 



+ e ab Ti N* H a H b 
+ M"ij k WWW. 



(1- 



We use the standard notation, which is explained in the footnote below. 3 The upper 
two blocks in Eq. (jl.8j) are R p conserving, the lower blocks violate R p ; the left blocks 

2 In 5*0(10) grand unified theories, the GUT symmetry can be broken by a Higgs field in the spino- 
rial 16-representation or in the 126-representation. If we exclusively use the latter, then i?-parity 
is conserved. If we include the former, then i?-parity will be violated after SO(10) breaking, see 
Refs. 1341 135] . From the low-energy point of view, the choice of breaking scheme, however, still 
appears arbitrary. 

3 H, Q, L represent the left-chiral SU (2) ^-doublet superfields of the Higgses, the quarks and lcptons; 
U, D, E, N represent the right-chiral superfields of the it-type quarks, d-type quarks, electron-type and 
neutrino- type leptons, respectively; an overbar denotes charge conjugation; a, 6, c and x, y, z are SU(2)w~ 
and S'?7(3)c-indices, respectively, k, I are generational indices, summation over all repeated indices 
is implied; S xy is the Kronecker symbol, e "' symbolizes any tensor that is totally antisymmetric with 
respect to the exchange of any two indices, with e 12 -- = 1. All other symbols are coupling constants, a 
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do not contain right-handed neutrinos, the right blocks do (hence they are absent in the 
models of Refs. [1-9]). Thus the superpotential of the MSSM is contained in the upper 
left block. We also consider non-renormalizable F-term operators of dimensionality five; 
having dropped all gauge and flavour indices they read 

LH U LH U , QQUD, QLUE, NLLE, NQLD, NUDD, 
QQQL, H V H U H V H U , UUDE, WNH V H U , NNNN, 

LH U H V H U , QH V UE, QQQH V , NQH U U, NQH V D, NLH V E, 

WNLH U , (1.9) 

grouped according to the same scheme as in Eq. (jl.8j) . Note that there are operators which 
conserve R p but violate baryon- and lepton-parity unlike the case with renormalizable 
operators. We shall focus on those operators that can lead to rapid proton decay, 

^QQQL, ^QQQH V - (1.10) 

M s is the string scale (~ 10 18 GeV), SI/' are arrays of dimensionless numbers. We do 
not consider the operators UUDE, because they contribute to proton decay only via the 
unknown quark-squark mixing matrix for the right-handed up-quarks, see for example 
Refs. 02IHJ|II|. 

This text is structured as follows. In Section El we review the idea of Froggatt and 
Nielsen within supersymmetry and consider some of its implications. In Section El we 
discuss the phenomenological constraints on the charges of the U(l)x symmetry. Then 
we parameterize and list the most general set of X-charges, consistent with anomaly 
cancellation and Eqs. (|l.llll.7jl . In Section we explain the procedure with which we get 
from the string scale where the model is naturally formulated down to the weak scale, 
as well as which approximations we utilize. In Section 5 together with Appendix we 
present and discuss our results. In Section 6 we offer a short summary. In Appendix [X] 
the models of Refs. [1-16] are discussed and listed, as well as the X-charge assignments 
of the other models. In Appendix El we demonstrate the validity of the mass matrices we 
use in this paper. 

2 The Scenario of Froggatt and Nielsen 

In the models we consider, the superpotential and the Kahlerpotential are expected to 
originate from a ten-dimensional heterotic superstring, see e.g. Ref. [13]: After compact- 
ification to four dimensions and below the string cut-off M s such a theory may exhibit 

priori they are arrays of arbitrary complex numbers (except Ay fe/A"^ being antisymmetric with respect 
to the exchange of the first two/last two indices, I\y and A"'ijk being totally symmetric with respect to 
any exchange of indices) . 
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at most 4 one anomalous local U(l)x (see e.g. Refs. |H1 EE] and references therein). We 
assume this anomalous local U(l)x indeed exists and has generation-dependent charges. 
For example, Xgi, i = 1,2,3, denote the U(l)x charges of the quark doublet superfields 
Q l and correspondingly for the other superfields. In addition to U(l)x, we assume an X- 
charged SM-singlet left-chiral superfield A, the so-called flavon superfield. For simplicity 
we choose the convention that Xa — — 1. 

Close to the string scale where U(l)x is unbroken, we have to replace the superpo- 
tential by a U(l)x gauge invariant extension, e.g. instead of G^ u \j Q l H u W we have the 
non-renormalizable operators 

g % q h u w I — J e [x Qi + x m + Xjjj] . (2.11) 

The complex couplings g^ u \j are assumed to be of 0(1), i.e. 

1 
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<\g { %\ <Vw. (2.12) 



Furthermore Q\x] is defined as 



1 for x being a non-negative integer , 

e[.r}={ ' ' ' (2.13) 

else (supersymmetric zero) . 

6[...] appears in Eq. (j2.11|) due to: (%) Holomorphicity of the superpotential: negative 
exponents of the superfield A are not allowed. The resulting texture is denoted a super- 
symmetric zero [49J. (ii) Fractional exponents are also forbidden, somewhat imprecisely 
the corresponding zeros shall be referred to as supersymmetric zeros as well. If an expo- 
nent of -4- is fractional and/or negative it shall be referred to as a naive exponent.^ The 
superpotential terms fiH v H u and KiL l H u are treated in exactly the same manner: 



. X H T>+ X H U 



X H T> + X H U 



■ X Li+ X HU 



Kj V H u — 9 



Xjj + Xjju 



(2.14) 



is a dimensionful parameter, the Ki are dimensionless and of 0(1). Note that the 
same mass scale is assumed in both expressions above. Likewise for the superpotential 
terms Hj N l and N l NK We again assume a common mass scale, in general different 

4 If there IS more than one anomalous U{1), one can re-express the charges such that only one ?7(1) 
is anomalous. For an explicit example see e.g. Ref. and references therein. 

5 However, as stated in Ref. [ST] one could "envisage fractional powers of the [super]field [A], stemming 
from non-perturbative effects" . 
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from MM: 



x- 



A*™ ll3 N l W ( A 



e 



x 



• JV* JVJ 



e 



NJ 



(2.15) 



is a dimensionful parameter presumably of C(M GUT ) or C(M S ) , the 7^ are 
dimensionless and of (9(1). If one is not willing to introduce yet another mass scale M^ N ' 
in addition to M 3 / 2 and M s (and possibly M^), one can set = M s and sufficiently 

suppress the righthanded neutrino masses by powers of e. We will explicitly deal with this 
point in Ref. (SOI- For a demonstration see the remarks about Ref. [10] in Appendix A, 
in particular Eq. (jA.ljl . 

Just like the superpotential, the Kahlerpotential also has to be replaced, e.g. instead 
of H(% Q 1 Q j we have 




e 



Xr 




X Q i X Qi 



e 



Xr 



x, 



QJ 



(2.16) 



The complex couplings h^\j are of 0(1) and form a positive-definite Hermitian matrix. 
For simplicity the gauge connection matrix is not written out as it is of no importance 
here. Note that these terms are not of the canonical form of the kinetic energy, which is 
important later on (c.f. Step 2 in Section 4). 

In addition to the expression above, we also obtain Kahlerpotential terms of the form 
— \ -x— 



H " ~M 



e 



Y 







* (if 9 


y 





— \-X rT ,-X hU 



a 



-V H' 



u 



a 



QiH u Ui 



M' 




e 



XflT) — XftU 



+ 



(2.17) 







X T 



-X qi ~X hU ~X- 



X H u 



e 



+ 



An» — X H u — X- 



+ 



together with the corresponding complex conjugate terms. The C_ symbolize the dimen- 
sionless 0{1) coupling constants. Z is an A-uncharged left-chiral superfield of the hidden 
sector, its F-term breaks supersymmetry. 

Just as with the U(l)x itself, the breaking of the U(l)x is also embedded in string 
theory. This can be traced back to the necessity to get rid of the mixed chiral anomalies 
of U(l)x with the SM gauge group. They are canceled by the four-dimensional field 
theory remnant of the Green-Schwarz mechanism j2H|, relying on two assumptions: 
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1. The left-chiral scalar complex dilaton superfield S shifts under a U (l)x gauge trans- 
formation, namely 

S - S + ^5 GS A x , (2.18) 

while the other superfields transform the usual way, e.g. H u — > exp(iX H u Ax) H u , 
Vx — > Vx + \ (Ax — Ax). Here, Vx is the U(l)x gauge vector superfield, Ax is 
the left-chiral superfield parameterizing the extended U(l)x gauge transformations 
and 5qs is a real number. 

2. The coefficient of the gravity-gravity-[/(l)x anomaly, ^4ggx, as well as the coeffi- 
cients for the mixed gauge anomalies SU (3)c-SU(3)c-U(l)x, SU (2) W -SU (2) W -U(l)x, 
U(1)y-U(1)y-U(1) x and U(l) x -U(l) x -U(l) x should fulfill, see e.g. Ref. 



-4ccx ! -4wwx i ^4yyx ! -4xxx ! ^4ggx ! 



2tt 2 ,5 G s. (2.19) 



kc kw ky 3 kx 12 

The fc... are the affine or Kac-Moody levels of the corresponding symmetry. The 
most easily constructed string models and thus so far most string models to date 
have 6 

kc = kw, kY = 5/3kw (2.20) 
The anomaly coefficients are given in terms of the X-charges as 

3 

^ccx = E ( 2 X Q l + X W+ X w) » ( 2 - 21 ) 



i=l 



Avwx = X H u + Xjjv + ^2 (3 Xgi + X L ^j, (2.22) 

1=1 

1 3 

•Ayyx = X m +X h t, + -J2(x q ,+8X w +2X w + 3 X Li + 6X-^j . (2.23) 

i=l 

Due to the possible existence of X-charged SM-singlets other than A and N l , 
the anomaly coefficients Axxx = TrX 3 and ^4ggx = TrX are not helpful to find 
constraints on the X-charges of the SM, so they are not listed. 

These two assumptions together ensure that the shift of the dilaton eliminates the anoma- 
lies of the X-charged current, because the dilaton couples universally to the topological 
terms e^ va ^ k,„ F~" F^h of the various gauge groups. For more details see e.g. Ref. 

6 The k... for non-Abelian groups have to be positive integers. For gauge theories arising from string 
theories fc... <?... 2 = 2g s , furthermore heterotic string theory always — no matter whether one has a 
GUT or not — gives rise to unification of the coupling constants, see Ref. jS2j. Hence at high energies 
gy = \/3/5 • gwi and thus ky = 5/3 kw- 
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To be entirely anomaly- free, the mixed anomaly U(l)y-U(l)x-U(l)x must vanish on 
its own, 7 

3 

-4yxx = X H u 2 — X h t> 2 + ^2 {xQi 2 — 2 Xjj- 2 + X-p- 2 — X L i 2 + X-g- 2 J = . (2.24) 

i=l 

It is crucial that Assumption 1 makes it necessary to have a modified Kahlerpotential for 
S, in order to have U(l)x gauge invariance: instead of being a function of S + S it has 
to be a function of 

S + S - 6 GS V x . (2.25) 

This "new" Kahlerpotential inevitably leads to string radiative corrections generating a 
finite U(l)x Fayet-Iliopoulos term G2E E3 EH] , its coefficient reads 

i FI x =l|^^ 2 ^GX. (2.26) 

£ M x tree level has to vanish in local supersymmetry, see Ref. Despite C, F1 x ^ we 

demand that (i) SU '(3)cxS77 ~(2) w xU (l)y and (ii) supersymmetry both remain unbroken 
at this scale, (i) is the case if the scalar components of all SM-superfields have vanishing 
VEVs. Taking this into account, (ii) is given if at least one of the X-charges of the 
remaining superfields whose scalar component gets a VEV has the opposite sign of £ FI x . 
This one superfield is the flavon superfield A, its scalar component acquiring the VEV 
(thereby breaking U(l)x) denoted by v oc y/^^x-i see a l so footnote ^1 Thus the string 
radiative correction is essential for the breaking of U(l)x, hence occurring near the string 
scale. This is why v/M s oc a/^4ggx/1927t 2 may not be a very small number. 8 On the 
other hand it must be ensured that v/M s is small enough so that we can identify 

-^- = e = 0.22. (2.27) 

Thus the V-charges should be moderately valued, otherwise ^4ggx would be so large 
(assuming that there is no fine-tuned cancellation among the V-charges) that it cannot 
be adequately suppressed by 1927T 2 . 9 

After the breaking of U(l) x we get from Eqs. (jl.8|) and (j2.11|) that 



X Ql + X H u + Xjjj . (2.28) 



Similarly, Eq. (|2.1fij) leads to 

H iQ \j = h {Q \j e lx Q*~ x Q> 1 e ia;„ - X, 



(2.29) 



7 Note that if the SM gauge group derives from an S'/7(5)-GUT, i.e. Xgi — Xjp — -X^r, X^i = X-jjr, 
and furthermore X h t> = —X H u 7 one automatically has Accx = Awwx = |-^yyx and Ayxx = 0. 

8 As can be seen from Eq. H2.26J1 . if U(l)x is non-anomalous (thus Aggx = 0) the string radiative 
correction vanishes; furthermore if in this case one chooses £ x appropriately one can break 

U(l)x at a much lower scale; such a model is sketched in Ref. 

9 As was shown in Ref. |3] (see also Ref. with Eqs. Ij2.19l rOSjl and using fc... g.. 2 — g s 2 , e is 

indeed of the correct order of magnitude. See also Ref. [50] , 
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and likewise for the other coupling constants. Depending on the U(l)x charges, e.g. the 
sum Xqi + X H u + Xjjj , we get generation- dependent exponents of e and thus generation- 
dependent suppressions of the G^ u \j . We thus have generation-dependent hierarchical 
coupling constants. It should be emphasized that the idea of Froggatt and Nielsen does 
not reduce the number of parameters, but only explains their hierarchy. As e is rather 
large, instead of the constraint in Eq. (j2.12j) we shall restrict ourselves to 



< 



(2.30) 



Due to the Giudice-Masiero mechanism, Kahlerpotential terms as in Eq. (|2.17J) also 
contribute to the effective superpotential, 10 giving 



M Mi 



3/2 



C H T> H U e lX HV+ X Hu\ 9 
C L i H U (:\ X L*+ X HU\ 



ix 



Mi 



3/2 



M-. 



3/2 



\XjjT> + X H U 

Xjj + Xfju 



(2.31) 



M: 



3/2 



a 



A" A" 



\X—+X—\ o 



IX- 



N< 



WW, 



M 



3/2 



C . A x qi +x hu+ x —\ 



l^Qi + X H u + Xjjj\ 



It is safe to ignore the contributions to the trilinear superpotential terms, as they are 
strongly suppressed by at least a factor of (^(lCr 15 ), 11 so that coupling constants as in 
Eq. ()2.28|) remain unchanged. This is not the case for the linear and bilinear terms: using 
the expressions above and Eqs. f)2. 14|. I2.15|) we obtain 

,— ,\ 2 



X 



A" 



+ M 3/2 M s C w e lx ^Q 



X- 



(2.32) 



+ M 3/2 Cj^ e lx TF +x w\ 



X^yj + X N j 



Ni 



(2.33) 



X^o + X^w 



+ M 3/2 C^fl" e l^+^l ©[(X^+X^ 



(2.34) 



Ki 



x L » +x 



+M 3/2 e l^+^"l0 



\Xj J i + Xj^-W 



(2.35) 



10 Supersymmetry is broken by the VEV of Fz, where Fz is the auxiliary field of Z. Thus one gets 
(Fz) 5 2 {6) = M3/2 M s S 2 (6). The ^-function turns the original Kahlerpotential term into a superpoten- 
tial term. 

11 See however Refs. They use the Giudice-Masiero mechanism to get small neutrino masses. 
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To summarize this section, the supergravity-embedded models which we consider 
are determined by the superpotential given in Eqs. (jl.8|) and (jl.lOj) together with the 
Kahlerpotential. The two potentials are generated close to the string scale by the scalar 
component of the flavon superfield acquiring a VEV. After U(l)x symmetry breaking, 
non-renormalizable operators generate renormalizable operators which are suppressed 
by powers of e times an unknown 0(1) coupling constant. The power of e is fixed by 
the generation-dependent quantum numbers of the superfields under the U(l)x family 
symmetry. One also has to take into account that via the Giudice-Masiero mechanism [4^ 
terms of the Kahlerpotential can contribute to the superpotential after U(l)x symmetry 
breaking. Numerically the effects are only relevant for the linear and bilinear terms. 



3 Phenomenological Aspects 

We would next like to derive a conclusion from the previous section (see also Ref. [65J). 
Consider Eq. (|2.19j) . it leads to 



~ (*4yyx + Avwx - | .Accx) = — Hp (k Y + kfa - ~ k^j ■ (3.36) 

Due to Eq. (l2~27IJ) the RHS equals zero, and with Eqs. (I2~2T1) . (l2~2^|) and (l2~2l?|) we obtain 
the following condition on the X-charges 

3 

X H u + X H v + ^ (x L i + — X Ql — Xjy[j = . (3.37) 

i=l 

Now the low energy mass matrix of the <i-quarks, M^ D \ is {H v ) Y^ D \ and the low energy 
Yukawa coupling matrix is the renormalization group evolved G^ D \ Disobeying 

virtually all of the steps which are outlined in the next section (to get from the string 
scale to the weak scale) we approximate here 

detY (D) = detg (D) ■ e 3X ^+^(^+^ , (3.38) 

with g \j = g( D \j @[Aq; + X h t> + Xjj-]. 12 The same holds for the mass matrices of the 
charged leptons. So from Eq. ()3.37j) we obtain 



detM(°) detgW 



(3.39) 



detM( £ ) detfl^) 
Assuming det g( D ^ ~ det g( E > and using 

= detM {E \ (3.40) 



12 Thus in order to have md ^ one needs detg( D ^ ^= 0. Using e.g. the simplest Ansatz for the g 



(D) 
ij ' 

namely all of them being equal to unity, is not an option. In Ref. [HDJ det g^ u ^ — was used to get a 
massless u-quark, see also footnote^. 
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Eq. (|3.39|) can be rewritten as 



e x HU+ x HT , = m d -m s -m b ^ 



With Eqs. (fTTj) . (1.2), and (1.3) we get 

X H u + X H v = or - 1 , (3.42) 

a result that is reproduced later on in a somewhat different way, see Table 1. Looking 
at Eq. (1231) we find that for X H u + X h t> = we have to have |M^| < M s , thus 
the hierarchy problem remains. On the other hand if X H u + X h t> = — 1 one has that 
H = e ■ M 3 / 2 . Hence M^' (which is also the mass scale for the Ki) is then naturally 
chosen to be M s . As explained in more detail below, the Ki can be rotated away which 
however changes the /i (c.f. Step 3 of the next section). So if ~ M s , one thus needs 
X L i +X H u to be either > 24 (so that ~ M s is adequately suppressed) or < and/or 
fractional so that the entry is forbidden. In that way the Ki do not generate an effective 
\x which is too large, see Eq. (|2.35j) . 

A similar kind of calculation was presented in Refs. [2H 122] : The conditions for 
Green-Schwarz anomaly cancellation, Eq. ()2.19|) . give 

^ *4yyx = - -4yyx + 2 -^wwx — -4ccx • (3.43) 

Expressing the RHS in terms of A-charges gives 

Y ^ 112 

- -4yyx = X h t> + X H u + {X^i + Xgi — + ~ 3"^^) 

i 

= 3X H v + Y,(X L >+X w ) + i (3 X H u + Y.( X Q* + X w) ) 

i % 

(3X H v+J2(X Q *+Xoi) )■ (3-44) 

i 

Similarly to the previous calculation, we get 

m e ■ m M • m r ( m u ■ m c ■ m t 
^ = < gP > 3 V <*">' ' . (3.45) 



2 
3 



m d -m s - m b 
{H v f 



Using Eqs. (1.1-1.5) and cot (3 = (H v ) / (H u ) this reads 

%yx m t 2 e 6 ^ ' 1 ' 2 or 3 )+(° or ^ , A N 

e s = ^ . (3.46) 

With (H u ) 2 + (F c ) 2 = (246 GeV) 2 one arrives at 

^Y_QC 6+(0.1,2 or 3)+(0 or 1) 1 + tal)2 P ( m t V /o 47 A 

tan 2 /3 v246 GeW ' 1 j 



11 



Experimentally, the parameters are restricted by [HEl EZ| 



2.4 < tan/3 < 50 and m t = (174.3 ± 5.1) GeV, (3.48) 

and therefore .Ayyx/5 ~ 6 + (0, 1, 2 or 3) + (0 or 1). It follows that .4yyx = is impos- 
sible, and hence (c.f. Eq. (J2.19)) ) ^4wwx = -Aqcx = are also not possible. So U(l)x has 
to be anomalous, and we need the Green-Schwarz mechanism, which fixes the breaking 
of U(l)x to be close to the string scale. 

Further phenomenological constraints on the X-charges arise more systematically 
from requiring the quark mass matrices to reproduce Eqs. (1.3), 1)1.4)) . 1)1.50 . ()1.6)) and 
(11.7)1 at high energies. There are four possible pairs of and (in Appendix IB1 

the validity of these mass matrices is illustrated): 

/ e 8 e 5 e 3 \ / e 4 e 3 e 3 \ 

G™ oc e 7 e 4 e 2 , oc e 3 e 2 e 2 , (3.49) 

V e 5 e 2 1 / \ e 1 1 / 

/ e 8 e 5 e 3 \ ( e 4 e 3 e 3 \ 

G™ oc e 13 e 4 £ 2 G (D) K £ 9 £ 2 £ 2 (350) 

V e 11 e 2 1 / V e 7 1 1 / 

/ e 8 e 6 e 4 \ / e 4 e 4 e 4 \ 

G (c/) oc e 6 e 4 e 2 , G (r>) oc e 2 e 2 e 2 , (3.51) 

V e 4 e 2 1 J \ 1 1 1 J 

( e 8 e 6 e 4 \ / e 4 e 4 e 4 \ 

G (C7) oc e 14 e 4 e 2 , oc e 10 e 2 e 2 . (3.52) 

V e 12 e 2 1 / V e 8 1 1 / 

The matrices in Eqs. ()3.49) 13.50)1 were suggested in Refs. [H Ej and more rigorously 
derived in Ref. [3 - 13 They are in accord with the lefthand choice of Eq. (1.7). Ref. [OH] 
pointed out that the righthand choice of Eq. (jl.7j) can be achieved from the mass matrices 
in Eq. ()3.51)) . 14 For completeness' sake, in analogy to Eq. ()3.50)) we found here that the 
matrices in Eq. ()3.52)) are also in accord with the righthand choice of Eq. (1.7). 15 

It is convenient to parameterize the X-charges in terms of the entries of the matrices 

13 Utilizing Eq. (|4.11(l . the matrices in Eq. (|3.5Q(I can be generated from naive exponents Xqi +X H u + 
Xjjj and Xqi + X h t> + X-^j given as follows, respectively: 



8 


-1 






-3 




13 


4 




• 


2 


1 


11 


2 


/,, 







/ 



The negative exponents give textures which are filled up in the process of canonicalization of the 
Kahlerpotential, see Step 2 of the next section. 

14 Note that G (c/) in Eq. lpT5l)l is the only matrix that is compatible with an SU(5)-GVT: G (u) 
is symmetric, being in accord with the 5't/(5)-requirement Xjp- — Xqi . To have furthermore X L i — 
X-jji, X-^; = Xqi one needs X L i = X L i = X L 3, y = 1, z = 0, in the notation below. 

15 Utilizing Eq. (|4.11(l . the matrices in Eq. H3.52(l can be generated from naive exponents Xni +X H u + 
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in Eqs. ll3.4HI3.52p : 



5 + y 3 + y 

X Qi +X H u+X w = I 7-y 4 2 | , (3.53) 

5-y 2 /;; 

4 + x 3 + y + x 3 + y + x 
X Qi + X H v + X w = \ 3-y+x 2 + x 2 + x \ , (3.54) 

1 — y + x x x 



with 



2 = 0,1,2,3 if y = -7,0,l and 2 = 0,1,2 if y=-6. (3.55) 

A negative y causes some naive exponents in , to be negative as well (see 

footnotes [H"] and IT""") and is thus forbidden; these textures are then filled up in the 
process of canonicalizing the Kahlerpotential. In addition to Eqs. ()3. 531 13*31)) . Eq. 1)1.1)) 
motivates the parameterization 

/ 4 + z + x \ 

X Li + X h t> + X w = 2 + x , (3.56) 

with ^ = 0, 1. Finally, Eq. (1.2) leads to 

X L a + X^j = Xqs + . (3.57) 

Solving Eqs. 1)3.53113.57)) and the conditions for Green-Schwarz anomaly cancellation, i.e. 
Eqs. (12.191 12 .21112. 24")) for the X-charges leads to the expressions given in Table 1; the 
free parameters chosen here are x, y, z, Xjji, X^2, X^. 16,17 In the Table in Appendix A we 
Xjjj and Xqi + X h t> + Xjjj- given as follows, respectively: 





As in footnote ^3 the negative exponents give textures which are filled up in the process of canonical- 
ization of the Kahlerpotential, see Step 2 of the next section. 

16 The reasons we have chosen to work with X L i,X L 2,X L s are as follows: (i) The quark sector is 
fairly well known, so Froggatt-Nielsen model-building does not need an X-charge of a quark to be an 
adjustable parameters, (ii) If G^ E \G^ N ^ and Y are without supersymmetric zeroes, one finds that the 
Maki-Nagakawa-Sakata matrix is approximately given by 

UMNS.^^-X^^ jusUike uC KM.. „ e \X Ql -X Qj \. 

furthermore (see Ref. [3 EI) 

17 If one believes that the actual expansion parameter is e 2 rather than e, one has to work with x = 0, 2, 
y = 1,-7, z = Q. 
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X R T> 






+ X L a (6 + 2x + 2z) - x(36 + 6a; + 5z) - 18y - 2^ 2 ] 


X H u 


= —XfjT> — 2 


Xqi 


= |[30- (X L i +X L2 + X L3 ) +3x + 6|/ + 4z] 




= Xgi-l-J/ 




= X Q i - 3 - y 


Xjjt 


= X h t> — Xqi + 8 + z 


Xjp 


= Xjjr-3 + y 


Xjp 


= X w -5+y 


x w 


= -X H v - X Q i + 4 + X 


X^2 


= X^r-l + y 




= X w -l + y 


X W 


= -X H v + 4 - X L i + x + z 


X & 


= -X h t> + 2-X L2 +x 


x w 


= —X H T> - X L 3 + X 



Table 1: The constrained X-charges. 
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have translated all the models given in Refs. [1-16] into the (standard) notation above. 
This should allow for an easier comparison. Note that one has 



-4ccx = 3 (6 + x + z) 
as was actually already apparent in Eq. (j3.46jl . and 



(3.58) 



Xjj. + XfjTl + X- 



4 + x + z 2 + x x 
4 + x + z 2 + x x 
4 + x + z 2 + x x 



+ X L i-X 



LA ■ 



(3.59) 



As an example, consider the first model of Ref. [2]. The charges are determined by 
x = 2, y — 0, z — 0, Xjj = —12, Xjj = —13, Xjj = 55. Thus one has, using Table 1 



XfjT> = 0, Xfju = 



and the quark and lepton charges 



Generation i 


Xq 1 




X Ui 






1 


2 
3 


16 
3 


22 
3 


-12 


18 


2 


1 
3 


13 
3 


13 
3 


-13 


17 


3 


7 
3 


13 
3 


7 
3 


55 


-53 



This results in the quark mass matrices which are given in Eq. ([3.49)1 . and (taking 
into account supersymmetric zeros) the leptonic mass matrix 

/ e 6 e 5 \ 
G {E) ~ e 5 e 4 . (3.60) 
V e 73 e 72 e 2 / 

4 The Procedure 

Taking a top-down approach we start off at or close to the string scale with a family 
symmetry, i.e. with a certain X-charge assignment for the left-chiral superfields. We 
determine the exponents of e of all allowed interactions, including the baryon- and lepton- 
number violating interactions. We thus also determine the supersymmetric zeroes. We 
then translate this into a low-energy model in order to compare with the measured 
fermion masses, their mixings and especially the experimental bounds on the barity- and 
lepton-parity violating coupling constants. To do so, we perform the following steps: 
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1. In models that contain right-handed neutrinos, the corresponding scalar superpart- 
ners can acquire VEVs which are determined by the minima of the scalar potential. 
The scalar potential is non-negative (apart from supergravity effects) and super- 
symmetry has not been broken yet (since we are just below the U(l)x breaking 
scale, thus quite close to the string scale), so the absolute minima of the scalar 
potential are given by the roots of the scalar potential. These roots in turn are 
given by the roots of the superpotential. 18 Finding the roots of W is equivalent to 
finding the linear constant shifts 

W(x,9) — >W(x,9)+A\ (4.1) 

such that the tree-level tadpole term Sj N l in the superpotential vanishes 19 [note 
that the Giudice-Masiero mechanism might produce, Hj ~ M 3 / 2 M S , c.f. Eq. (|2.32|) ]. 
This shift results in constant and thus harmless terms (ignoring the problem of the 
cosmological constant) in the Kahler- and superpotential, as well as the redefinitions 
of the superpotential coupling constants: 

K t — Ki + A j , 

H — > fj, + T, A\ 

Si — > Si + 2 r y A 1 + 3 A"' ijk A j A k = 0, 

r y — > Tij + 3 M" ijk A k . (4.2) 

It is important to notice that the equations above lead to a mixing among the 
coupling constants. 

2. The breaking of U(l)x generates a Kahlerpotential which does not have the canon- 
ical form jlH] ■ Thus we must perform a transformation of the relevant superfields 
to the canonical basis [7T]. For example, for the quark doublets we obtain for 
the relevant Kahlerpotential term (where the H^\j are Kahlerpotential coupling 
constants which are to be canonicalized) 



Q l H ( % Q 3 = \^/D HiQ) U H(Q) Qj [y/D HlQ) U H(Q) Qj . (4.3) 

D H (q) is a diagonal matrix, its entries are the eigenvalues of the Hermitian matrix 
ffC); the unitary matrix U H (q) performs the diagonalization. We define the matrix 



C iQ) = ^/D^U HW . (4.4) 

18 Of course this consideration neglects the part of the scalar potential which arises from the £)-terms. 
However, this is justified by our starting point: We are working with one flavon field, i.e. U(l)x gets 
broken predominantly by the scalar component of A acquiring a VEV, rather than the scalar components 
of the iV* acquiring VEVs. This is equivalent to u > A 1 , which is justified with hindsight by the A* 
being functions of 3, , . 

19 Si is not regenerated by the RG-flow after having been shifted away, see the expression in Ref. |7(1] 
for the beta-function of a general linear term of a superpotential. 



16 



The transformation to the canonical basis can then be written as 



Q l -> C {Q) \ Q j . 



(4.5) 



The redefinition also affects the superpotential, e.g. the u-type quark Yukawa 
couplings 

G (U) _^ 1 C (Q)- lT q(U) c (U)-\ (4 6) 



In this way the canonicalization of the Kahlerpotential "fills up" possible super- 
symmetric zeroes (which are due to holomorphy and fractional X-charges [4"9]). 20 
in this particular case in the mass matrix 

Since L\ H v have the same gauge quantum numbers under SU(3)c x SU(2)w x 
U(1)y, they can be mixed by kinetic terms. So the canonicalization of the Kahler- 
potential causes mixing of G^ E \j with A^y, G^ D \j with A'jjy (an example is given 
at the end of this section), /i with Ki, G^ N \j with Tj, and ^ijki with ty'ijk- 

To maintain the canonicalized form, all further field redefinitions have to be uni- 
tary. 21 To demonstrate an exact canonicalization, consider the quark doublets with 
the U(l)x flavour charges Xqi = 4, Xq2 = 2, Xq-a = and setting h^\j = 1 (c.f. 
Eq. ()2.16Jl ). For the Kahlerpotential matrix of the Q*-superfields we obtain 



H (Q) 



e 4 e 2 



(4.7) 



This is canonicalized by the matrix (c.f. Eq. (|4.4)P 





/ 


i 





1 

V2 


C (Q)-t_ 1 




MJ4. 


w+f- 










2 


~2~ 




V 




w-f+ 
2 


Ul_ 

~2~ 



(4.8) 



with w± 



1 ± 



3^ 



Vs+^ ' /± = As e is a small quantity, the expres- 

sion above can be approximated to leading order in e as 



C (Q) 



-it 







i 

V2 



\ 



1 , 3 6 2 

2 ^ 8x/2 



V 



3_£ 

8^ 



1 

_L 



5 e 2 



I _|_ 3 e 2 
2^ 8V2 



5 e 2 
8 



(4.9) 



SV2 / 



20 There are exceptions though: If one works with certain X-charge assignments that involve fractional 
charges (e.g. the one stemming from Xj^ = j and X h t> = |) the four zeroes of are not filled up. 

This seems to have first been mentioned in Ref. @]. 

21 Actually the canonicalization is affected by loop-diagrams, but these effects can be neglected because 
they are small: e.g. the graphs contributing to the propagators of the neutrinos are similar to the graphs 
that radiatively generate tiny neutrino masses. 
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However, while this article was reviewed, a preprint [7.2] came out, stating and nicely 
demonstrating that in many cases one can utilize a further unitary transformation 
to bring 1 to its "standard" form, see Eq. ()4.11|) . 

3. Next we perform a unitary field redefinition of the fields (H v , L l ) such that the term 
KiUH u vanishes. This was first realized in Ref . [Z2] , however, the rotation matrix 
was approximate and in addition only valid for real fi and Ki, with Ki/ fi <C 1. 
We have implemented the full transformation for arbitrary complex K^fi. With 
K = \/Ki*K l and Ai = a/ + K 2 (summation over repeated indices implied) 
the exact transformation is given by 



/ 



M 



K 




(4.10) 



The SU{A) matrix acting on the fields (H v , L 1 ) is given by Ufc^ (we shall refer to 
it as a rotation matrix, since for real coupling constants it is an SO{4) matrix). The 
rotation gives further mixing of the coupling constants: G^ E \j <-» A^y, G^ D \j <-> 
A' ik j, G^ N \j <-> Tj and ^ijki <-> ^'ijfc- As explained in Section if rrie/m^ fa e 5 , 
one naturally has ~ M.. So one obtains the phenomenologically unacceptable 
relation that fa-M is much larger than the weak scale, unless the Ki are suppressed 
by at least e x v +x H u ~ e 24 (e 24 M s ~ weak scale), see the text below Eq. (J312J) . 

4. In the next step, the coupling constants are evolved according to their renormal- 
ization group equations (RGEs) down from the £/(l)x-breaking scale to the weak 
scale. The one-loop RGEs of the ftp coupling constants are given in Ref. [Hj, the 
two-loop RGEs of the gauge coupling constants and the coupling constants of the 
superpotential of the fi p -MSSM were derived in Ref. [7J3> the one-loop RGEs of 
the soft supersymmetry-breaking coupling constants of the fi p -MSSM are given in 
Ref. [33], see also Ref. [TH] for the third generation only. The RG-flow once again 
causes G^ E \j <-> A ik j, G^ D \j <-> A' ik j and G^ N \j <-> Tj mixing and regenerates Ki 
[77]. Likewise for the corresponding soft coupling constants. 

5. The superfields U and H v are rotated unitarily such that only the scalar neutral 
CP-even components of H v and H u acquire VEVs. The relevant matrix is given 
by Eq. (|4.10|) . with the appropriate replacements. The rotation again mixes the 
regenerated K % with /i, and again causes G^ E \j <-» A ik j, G^ D \j <-» A' ik j, G^ N \j <->• Tj 
and ^ ijk i <-> ty' ijk mixing. 
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6. In order to compare with experiment, we must transform into the mass basis of the 
quarks, which once again skews the coupling constants. As we will assume massless 
neutrinos here, we do not have to be concerned about the mass basis of the leptons. 
The inclusion of neutrino masses will be discussed elsewhere. 22 

7. Finally, in the last step, the tree and loop contributions to the mass matrix of the 
neutral fermions are determined. The neutralinos, and the regenerated operators 
KiL l H u can cause non-zero neutrino masses. The ^t p operators allow loop graphs 
that radiatively generate masses for the neutrinos, see e.g. Refs. jZHl EUl IHD IH2]- 23 
In this paper we have not taken into account phenomenological constraints on the 
A-charges that arise from the neutrino sector. 

Note that if all of the (highly undesired) couplings A"^ are initially forbidden due to 
naive fractional and/or negative A-charges, none of the above steps can regenerate them. 

Of course, performing all these steps entirely satisfactory requires the precise knowl- 
edge of the unknown a priori arbitrary G(l) coupling constants g^ U \j, h^ U \j, C h t> h u, etc. 
Thus the only thing we can do is assume that there is no fine-tuning leading to accidental 
cancellations and to perform order-of-magnitude estimates. We would now like to point 
out several caveats and also comment on how our procedure differs from relevant previous 
work: 

1. Step 1 has been neglected in the past; we have taken it into account whenever 
A^pr, A-^2, Xjjj were given (this is not the case in the model presented in Ref. [10], 
the models presented in Ref. [12] and the model in Ref. [16]). The right-handed 
neutrino charges did however not change whether the model under consideration 
is compatible or not: the Xj^ in model [13] are fractional and thus no HjX* is 
generated, the models [14] and [15] are already ruled out even if the explicitly given 
Xjr? are not considered, likewise for the models [17]- [21]. 

2. For Step 2, in the past the approximation 

C«)-y~el*tf-*tf1 (4.11) 

was utilized (likewise for the other fields), see Ref. |2]- 24 For example, for the choice 
of quark charges Xqi =4, Xq2 = 2 and Aq3 = 0, the transformation matrix is given 

22 The neutrino mass matrix pattern required for the solar and atmospheric neutrino problems can also 
be obtained via the idea of Froggatt and Nielsen, see e.g. Ref. |78| . 

23 The corrections to the neutrino masses are totally negligible if the model contains right-handed 
neutrinos, in the same way as the corresponding corrections are unimportant for the masses of the 
charged fermions. 

24 This result can be reached as follows: Ignoring the fact that \Xq* — Xgj | could be fractional, 

see Eq. I|4.7|). Now \X Qi - X Qj \ < \X QJ - X Qk \ + \X Qk - X Qi \, from which follows that 
£ \ x q'- x qj\ > e \ x Qi- x Q k\ e \ X Qk- x Q i\ j an( j hence, neglecting prefactors of 0(1), ^ x q-~ x q^ w 
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by 

(I e 2 e 4 \ 
= e 2 1 e 2 . (4.12) 
V e 4 e 2 1 ; 

This is in gross disagreement with the exact result in Eq. ()4.9j) . Of course, Eq. (|4.fJj) 
relied on a fine-tuning (all h^\j=l), however it illustrates that Eq. (|4.11j) might 
be a very poor approximation. The correct matrix should be found somewhere be- 
tween these two extreme cases (see however the text below Eq. ()4.9|) and Ref. [HI]). 
Keeping in mind that the canonicalization of the Kahlerpotential is very important 
for filling up supersymmetric zeroes, one sees that Eq. ()4.9|) skews the coupling 
constants much more than Eq. 1)4.11)1 and thus leads to a much stronger unwanted 
attenuation of possible hierarchies. In this paper we shall utilize Eq. (j4.11j) . being 
the more conservative approach. 

3. Step 3 has been explicitly mentioned in the past only by Refs. [3IH|- We utilize an 
approximation to leading order in e. 

4. For the time being, we shall use the most up-to-date GUT-scale bounds on the 
single fi p coupling constants given in Ref. |HSj- For products of two fi p Yukawa 
coupling constants we use the weak scale bounds 25 of Refs. ISSHHHESHHEl (instead of 
using some of the bounds of Ref. [HI] one might wish to use the bounds of Ref. |87j : 
they are stricter, however less conservative). The bounds are compared in units 
of (m/100 GeV) 2 , m being the mass scale of the scalar particles with R p = — 1. 
For the bounds on higher- dimensional operators we followed Ref. [21 IS]- They 
state ^n2i < 10~ 8 m/100 GeV, deriving from the experimental bounds on proton 
decay. To take into account the large uncertainties in determining this bound (see 
e.g. Ref. jHH]), we worked with the more conservative 

^lm < 10~ 7 - — . (4.13) 

100 GeV V ; 



In addition, we always included a factor of tolerance 1 / ■s/e ~ 2 to take into account 
numerical prefactors C{1). We initially set m = 100 GeV for all bounds and then 



J2e\ X Q^ X Q k \ e |X « fe ~ X Q l1 . Thus 
k 

y^Qi e \ x Q*- x Qi\Qi & rqi e \ x Qi- x Qk\ e \x Qk -x Qi \ Qj ^ 

and hence, comparing with Eq. I|4.3|l . h^q^-^q^ C^\j. Furthermore, the el x <5« _x ejl form a matrix 
that is approximately equal to its square and thus approximately equal to its cube, so that this matrix 
is approximately equal to its own inverse, so that 



^(Q)- 1 ! M f \x Q *-x Q1 \ = e \x Qj -x Qi \ _ r^wr 1 l = r C (Q) 

-ij L Jij 



-it 



This approach is conservative, as the bounds get tighter when going to higher scales. 
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indicate whether a model is not ruled out if we allow 100 GeV< m < 500 GeV. 
Heavier sfermion masses lead to weaker bounds. 

5. Concerning Step 5, we shall work with the approximation that the VEVs of the left- 
handed sneutrinos vanish as has always been done in the past (except in Ref. 0). 
This is motivated from low-energy phenomenology: the coupling constant of the 
soft supersymmetry breaking bilinear term with a scalar L 1 and a scalar H u has to 
be very small compared to the Higgs-VEV, and the soft supersymmetry breaking 
4x4 mass matrix of the scalar (H v , L l ) has to be nearly diagonal. Hence the 
VEVs of the left-handed sneutrinos have to be very small, see Ref. [75] and also 
Ref. [35] . When considering the neutrino masses generated through the K$ in order 
to be consistent we must then include the VEVs for the scalar neutrinos. 

6. The matrices that biunitarily diagonalize the quark mass matrices, namely U^ Ul \ 
U^ Ur \ U^ Dl ^ and U^ Dr ' are given in Ref. jSH] as a function of the entries of the 
quark mass matrices, as well as JJ 01 ™ = xj( u l)^ij{Dl) _ We shall not state the 
corresponding expressions here, because they are long and complicated, being a 
function of the entries of the quark mass matrices. Since we do not know the 0(1) 
coefficients and are hence only performing order-of-magnitude estimates, we set all 
coefficients equal to 1; however, to avoid accidental cancellations we thus have to 
replace all minus-signs in the expressions in Ref. [EH] (this is also the procedure 
adopted in previous work). We have checked that the corresponding results are in 
good agreement with the exact results, having chosen a randomly generated set of 
(9(1) coefficients (see Appendix 2). 

7. Due to the uncertainties on masses and mixings in the neutrino sector we shall nei- 
ther use radiative corrections to neutrino masses (based on fi p coupling constants) 
nor Dirac- and Majorana-masses (when considering right-handed neutrinos) to con- 
strain the V-charges. 

It is important to note that all simplifications which we have adopted here lead to a 
more conservative approach than calculating exactly (assuming that there are no acci- 
dental cancellations). 

Now the various transformations (apart from renormalization group flow) in Steps 1-7 
are all unitary or at least almost unitary. Thus the determinant of the mass matrices 
should not change substantially, so that apart from a few accidental cancellations and 
"accidental adding-ups" (and apart from supersymmetric zeroes being filled up) the en- 
tries of the mass matrices (and also the entries of the other Yukawa coupling constants) on 
the average keep their order of magnitude. In order to estimate the possible "accidental 
adding-ups" in our conservative approach, we have calculated the effects of Steps 1-7 on 
the coupling constants to higher order in e (which can contribute quite a lot as e is not a 
very small number). We used the approximate expressions of the various transformations 
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as for example Eq. (j4.11j) but kept all powers of e in the result. Consider e.g. the charge 
assignment of Ref. [2] which is displayed at the end of the last section. G^ u ' is given by 
Eq. fl33H|), and C^) -1 , C&V are given by (c.f. Eq. (|4.11|) ) 



e 
1 

'2 



(4.14) 



G<& is then (c.f. Eq. 



9 e 8 
6 e 7 + 3 e 9 
3e 5 + 3 e 9 + 3 e 11 



6 e 5 + 3 e 



ii 



3 e 3 + 3 e' 



3 e 



13 



4e 4 + 2e 6 + 2e 10 + e 12 



2e 2 + e 4 + 2e 6 + e 8 + 2e 12 + e 14 



2e^ + 2e b + 3e» +e 12 +e i4 1+ 2e 4 + e e + e 8 + 2e iU + e i4 + e 



,14 



10 



,14 



.16 



In the past (and in our conservative approach as well) all numerical prefactors and all 
higher orders in e were neglected. If they are taken into account one gets with e ~ 0.22 
that 



G (U) 



£ 6.5 e 3.8 £ 2.3 
£ 5.8 £ 3.1 gl.5 
£ 4.3 e 1.5 1 



(4.15) 



differing drastically from the conservative result, which is the same as in Eq. ()3.49|) . This 
is not meant as a demonstration that the conservative approach produces the wrong mass 
matrix, it is just an indication that some of its entries potentially can get quite large due 
to "accidental adding-ups". 

Using the same X-charge assignment as before, we give an example which illustrates 
how important it is to take into account the U <-> H v mixing when canonicalizing the 
Kahlerpotential which has been ignored in the past. For the Ifi p coupling constants we 
obtain 



A' 



ilk 



0, A' t2k = 0, A', 



?>k 



.55 



(4.16) 



The model has no right-handed neutrinos so that we do not have to perform Step 1 which 
is outlined above. To canonicalize the Kahlerpotential (Step 2) the basis transformation 
for the H v and L l fields is given by 





/ 1 


e 12 


£ 13 


e 55 


\ 


1 


e 12 


1 


e 


£ 67 






£ 13 


e 


1 


.68 






V e 55 


£ 67 


e 68 


1 


/ 



(4.17) 



This matrix mixes G^°> with the A' given in Eq. (|4.1fij) : 









V e 55 GW J 



12 r 1 iD) i Y22q(D) 



/ 



e 13 Qip) + f 123g(D) 

\ 2 e 55 



(4.18) 
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Dropping higher order terms in e and neglecting 0(1) prefactors this gives 

/ \ 
e 12 GW 

e 13 Q{D) ■ 

V e 55 J 

With e = 0.22 we conclude that the largest A'ij k is of 0(iCr 8 ). This leads to rapid proton 
decay, because it turns out that the largest A"^ is of O(10~ 7 ). So the model has to be 
regarded as incompatible with broken .R-parity, which is why it is labeled as "no" in the 
Table in Appendix A. Had we not taken into account the L % <-> H v mixing (hence G^ D > 
not mixing with A'^ ) we would have obtained 

/ \ 

e 122 Q(D) 
e 123 Q(D) 

\ e 55 ) 

instead, giving A'^ < O(l0~ 38 ), which is in agreement with current experimental bounds 
on proton decay. The next step in the canonicalization of the Kahlerpotential would be 
to replace G (D) by C (Q)_1 -G iD) -C^ D) , and likewise A' ijk , but we shall not go into 
the details anymore. 



5 Results 

We have taken the existing Froggatt-Nielsen models given in Refs. [1-16] and embedded 
them into our general framework. As explicitly given in Table 1, after taking into account 
all the constraints from the quark masses and mixings, and the lepton masses, as well 
as the Green-Schwarz anomaly cancellation conditions, all U(l)x charges and thus all 
models can be expressed in terms of three integer parameters and three U(l)x charges 
which we choose to be (c.f. footnote tt^j) 

x, y, z G 7L, and X L i, X L 2, X L s . (5.1) 

The constraints on the integer parameters are given in Eq. (|3.55jl . as well as z = 0, 1. 
The specific choice fixes the quark mass matrices given in Eqs. (|3.49J) - (|3.52|) as well 
as the charged lepton mass matrix. As a first step, in columns 2-7 of the Table A in 
Appendix A we determine the parameters x, y, z and the charges X L i for all models in 
the literature. This makes it possible to easily compare them. In Appendix A, we give a 
detailed discussion of the models and how the charges were reconstructed out of the given 
information. In Table IXj we have models with small integer charges (< 10) as in model 
[1], the first model in [4], the models in [5], [7], [14], and [15]. Models with fractional 
charges constitute the rest. Note that X L i, X L 2, X L a being integer is only a necessary but 
not sufficient condition for all X-charges of MSSM particles to be integer. In order to see 
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this, it is important to go back to Table 1 and observe that only a special combination 
of x, y, z and the -charges guarantees that Xqi and X h t> are integer. This is why 
e.g. the third to last model in Ref. [3] does not have purely integer X-charges, although 
X L i = 6, X L 2 = 1, X L 3 = 4. 

Having defined a unified formalism for all models, we developed a numerical code 
which implements the steps outlined in Section 0] for an arbitrary model defined by the 
choice in Eq. (|5.1jl . We thus determine for each model the order of magnitude of the 
i?-parity violating Yukawa couplings and the dimension-five operators at the GUT scale. 
After running the renormalization group equations one then obtains the corresponding 
weak scale predictions. Of course this occurs under the assumption that there is no 
additional discrete symmetry such as i?-parity or baryon-parity which forbids a subset 
of terms. We can then compare these predictions with the existing experimental bounds 
on the .^-couplings and the dimension-five operators and thus either allow or exclude a 
given model of fermion masses. 

For single i?-parity violating couplings the weak scale bounds have been converted to 
GUT scale bounds in Ref. jHSl using the RGEs of Ref. [ZHj- We thus compare directly 
to these. However, we necessarily predict more than one non-zero i?-parity violating 
coupling at the weak- or the GUT-scale and we must take into account the bounds on 
products of couplings. These can often be much stricter than bounds on single couplings 
due to extra flavour changing neutral current effects [H31IE3IEI3IH2] and of course the strict 
bounds from proton decay [HE] • In column 8 of Table [X] in Appendix A ( "Compatible 
with Exp.") we list whether the model satisfies the experimental constraints on baryon- 
and lepton-number violation or not. 

• u yes" means that the model is compatible with the experimental constraints after 
performing the procedure outlined in Sect. H| 

• u (no)h. " means that the model is compatible if one calculates conservatively, but 
not compatible if one takes into account the summation of higher orders in e (see 
the example in the previous section). A moderate amount of fine-tuning such that 
higher order terms cancel would make the model compatible. 26 

• "no" means that the model is incompatible even to lowest order in e. 

• "yes/ no" means that the model is compatible if one uses the freedom to choose X^3 
correspondingly, but the model is not compatible if the particular value of X L z is 
used which is given in the table. 

• The models denoted by * change their status when the relevant sfermion mass is 
increased. The details are given in the table caption. 

26 Strictly speaking taking higher orders into account also rules out the mass matrices from which one 
starts out, see Eq. (|4.15|) 
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5.1 Discussion 



We find that only three models from Refs. [1-9] (the models without right-handed neutri- 
nos) pass with flying colors, i.e. produce a "yes", (i) One model of Ref. [2] is effectively 
i?-parity conserving due to the very high value of X L 3. (ii) The model in Ref. [S] is com- 
patible due to complicated fractions for the X-charges. (Hi) The model |5J also survives 
due to complicated fractional charges. 

Many models, like most of the models in Ref. [7], are incompatible solely due to the 
bounds on higher dimensional operators. As an example, consider the first model of 
Ref [7]: x = y = z = 0, X L i = X L s = —6, and X L 2 = —3 give Xqi = 5, Xq2 = 4, so 
that the operator Q 1 Q l Q 2 L 3 is suppressed by e 8 , which is a factor ~ 55 too large. 

It is important to note that including the exact transformation matrix to the canonical 
basis for the kinetic terms, C ", as well as the exact matrix for rotating away the LH U - 
term makes it even harder for models to survive. 

In the case without right-handed neutrinos, we see that a flavour dependent U(l)x 
cannot solve the problem of baryon- and lepton-number violation if we demand small 
integer or simple fractional X-charges. We have thus also considered the case of an 
additional discrete symmetry: i?-parity, baryon-parity, or lepton-parity |2Ej. In the last 
column of Table 1X1 we show which set of couplings must vanish for the model to survive. 
Thus for example the second set of models in Ref. [1] survive (for certain n) if \1/ = A" = 0, 
which can be guaranteed by baryon-parity. Similarly, the third model in Ref. [3] is 
compatible with the experimental constraints if A = A' = 0, which can be achieved by 
lepton-parity (anomalous). The first model in Ref. [7] becomes compatible if ^ = A' = 0, 
which is also obtained by lepton-parity. In this manner, we can obtain models with small 
integer charges, but we also lose a central part of our motivation. If there is no entry in 
the last column the model is either compatible or beyond repair. 

When we allow for right-handed neutrinos, Refs. [10-16], we find a significant number 
of allowed models. Here the X-charges of the right-handed neutrinos present extra degrees 
of freedom, making the charge assignments more flexible to avoid bounds on the fi p 
coupling constants. (If one does not have right-handed neutrinos one is forced to have 
non- vanishing A^ and/or A'^ in order to get neutrino masses.) However, again we 
find that the models with integer X-charges ([14] and [15]) are not compatible. The 
models [11,12] automatically have conserved i?-parity and no QQQL operator, i.e. \& = 
0, due to the fractional charges (although X l i,X l i,X L 3 are integers). On the other 
hand, while the models in [10,2] certainly can fulfill all constraints if one chooses an 
appropriate X^s, certain choices of X^3 providing special realistic cases are in conflict 
with the constraints. The charges are not "fractional enough" , that is, they still produce 
positive-integer exponents of e for certain $ p coupling constants, thus not forbidding ]jtp 
by fractional X-charges as is the case in Refs. [11,12]. 

For completeness in Table B we also investigated and listed the Froggatt-Nielsen 
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models (one in Ref. [16], the rest in Refs. [17-21]) which exhibit a slightly different 
fermionic mass spectrum than the models in Refs. [1-16]. 

6 Summary 

We have presented a systematic embedding of four- dimensional Froggatt-Nielsen models 
with a single Froggatt-Nielsen field into supersymmetry. We have developed a simple 
notation into which we have translated all models [fulfilling Eqs. (1.1-1.7) and the Green- 
Schwarz conditions] which we found in the literature. We have then extended (where 
necessary) the models to include their predictions for baryon- and lepton-number viola- 
tion, including the dimension- five operators. Throughout, we have taken into account 
the modification of the fermion mass matrices due to canonicalization of the Kahler po- 
tential, and the rotation of the H v and Lj fields. We also consider the case of summing 
sub-leading powers of e, which can lead to significant effects. We have outlined in detail 
how the models are translated from high-energy textures to low-energy predictions. We 
have discussed the required steps and pointed out where we have improved on the exist- 
ing literature. In our Table in Appendix A, we show how many of the existing models 
fair when compared to the most recent bounds on B- and L-parity violating coupling 
constants. Only very few models survive unscathed. In the future we hope to improve 
on our procedure as mentioned in the steps. We also wish to extend our discussion to 
neutrino masses which we have not yet taken into account. 
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A The Catalogue of Froggatt-Nielsen Models 

The first nine references deal with models with no right-handed neutrinos, so no 5< N l 
term has to be shifted away. The values for the six parameters (x, y, z, X]j) were ob- 
tained as follows (the references are listed chronologically according to their appearance 
on the arXive): 

• P explicitly states four sets of A-charges, three of these sets fit the class of model 
considered in this paper (in the third model there is a typo in the A-charge assign- 
ments: instead of ei = — 1 + q, = 5 + q it should read e\ = 5 + q, = — 1 + q). 
The construction of these models did not take into account the Giudice-Masiero 
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mechanism, higher-dimensional operators and ftp. It turns out that the models are 
in conflict with the bounds on \l/ independent of the value chosen for X L 3. For low 
values of |X L 3|, the models also disagree with bounds on Ijt p . 

|2j states two complete sets of X-charges. The authors take into account higher- 
dimensional and ^-operators as well as the Giudice-Masiero mechanism. The 
second model is in accord with all constraints due to the large value of \X L a |. The 
first model was treated as an example in Section 4. The paper presents a further 
model with right-handed neutrinos that is also discussed in Ref. jTUj, see below. 

jHj considers the choice y = —6, (in our notation). The corresponding quark mass 
matrices must be combined with four different explicitly stated charged lepton 
mass matrices. These are obtained through four different choices for the charge 
differences X^x — X^s and X^i — X^3 . The first model is given (in our notation) by 
x = 0, the second one has y = —6 with x = 0,1,2. The third model is required to 
have x = 0, and for the last one x = 0, and y = —6 is demanded. In addition, the 
authors work with £q = X^a + Xqs + X-^a to be 4 or 5, —4 or —3, 7 or 8, and —7 
or —6, respectively. They thus have 2 • (2 + 3 + 2 + 1) = 16 X-charge assignments. 
All are such that no dangerous higher- dimensional operators arise. But the authors 
do not take into account the Giudice-Masiero mechanism, the LH u -term, or the 
mixing of L l and H v . They do consider ft p (explicitly needed to generate neutrino 
masses), however twelve models disagree with the experimental bounds at lowest 
order in e and the four other models disagree when including higher order terms. 

[3] explicitly states the three sets of X-charges. The authors take into account 
higher- dimensional operators as well as U <-> H v mixing and thus mixing of QQQL 
and QQQH V . The first two models are in accord with the bounds to lowest order 
in e, however not so to higher order. The third model is ruled out because it does 
not agree with ft p : X L i + X L z + X E i = and thus Ai 3 i = 0(1). 

[2] explicitly states an incomplete but sufficient set of X-charges. They disregard 
higher- dimensional operators, and the corresponding constraints are not met. 

jH] explicitly states two sets of X-charge assignments which however are not compat- 
ible with the Green-Schwarz mechanism. But from their requirements (translated 
to our notation) of y — 0, z = 1, X L i — 2 = X L 2 = X L z and 7 < X L a — X h t> < 9 
and tan/3 ~ 50 =^ x = one can extract one X-charge assignment which is 
compatible with the bounds. 

[7j states 66 (!) X-charge assignments, many of them however referring to quark 
mass matrices that deviate slightly from the ones considered here: The naive expo- 



27 



nents in and are given as follows 





We shall not consider these charge assignments. Focusing on Eq. ([3.49)1 . we are 
left with 26 X-charge assignments. Note that their 2 nd , 12 th and 15 t?i model are 
special cases of the 1 st , 2 nd and 2 nd model of Ref. [Tj, respectively. Their 8 th model 
equals the one in Ref. [5] . The authors did not take into account higher-dimensional 
operators; the corresponding bounds rule out many models. Although the authors 
do consider ftp, some of the models are ruled out due to disagreement with more 
recent bounds on ft p coupling constants. In the end, none of the models survives 
even to lowest order. 

|Hj explicitly states a set of X-charges. All constraints are carefully obeyed, leading 
to a valid model. The price to be paid are highly fractional charges. 

jH] explicitly states a set of X-charges. The author does not take into account 
higher-dimensional operators, as the corresponding bounds rule out the model. 

[TUj and [2j (the first reference is a talk discussing one of the models - this time with 
right-handed neutrinos - of the second reference) demand X L i — 3 = X L 2 — 1 = X L 3. 
This still leaves X L 3 undetermined. Combining this with the fact that the authors 
work with y = —6, 0, we get with z — 0, 1 that there are 2 x (3 + 4) = 14 charge 
assignments with X L z being a free parameter. Although explicitly including right- 
handed neutrinos, no corresponding X-charges are stated. Of course one might 
argue that it is always possible to choose a value for X^3 such that .R-parity is 
exactly conserved or at least broken in a way that experimental limits are obeyed. 
However, to consider a special realistic case, the authors demand the Dirac and 
Majorana neutrino mass matrices to be without supersymmetric zeroes and the 
(3, 3)-component to dominate, i.e. one needs Xj^+Xjju+Xjj > X-^+X H u+X L z > 
and X^ + X^y > X w + X w > 0. From Eq. (14) of Ref. QU|, we thus find 

(H u ) 2 e 2{x * u + x ^ 
m VT ~ ^ L ^ • (A.l) 

With y = 0, —6 the X-charge assignments are not GUT compatible, one thus 
assumes M = M s ~ 10 18 GeV (and not M = M GUT ~ 10 16 GeV). The authors 
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demand m Vr = 0.1 eV. Combining this with e ~ 0.22, (H u ^ ~ 100 GeV, we get 
that X H u + X L 3 ~ —3; we can thus fix X L z. 

[TTj and (the first reference is a talk discussing one of the models of the second 
reference) present a model (aimed at atmospheric neutrinos) that is based on the 
quark sector in Ref. j2| (Ref. ^2] is a kind of continuation of Ref. j2j. Ref. JJJ 
also reviews the model of Refs. [21 CD]), so y — —6,0. The authors demand | < 
X L 2 = —X L 3 < |, X L i — X L 2 > 3, which they specify to Xf,i = 5,X £ 2 = 2, giving 
X^r— X^" = 2 — 1. As they furthermore suppose Xp- = X^-, the authors get z = 1. 
This leads to m e /m T ~ e 5 , as is correctly stated in Ref. [12] . It should however 
be mentioned that both Refs. JJJ and J2| explicitly demand X#« = X^-d = 0, 
which is in contradiction to 2 = 1. Furthermore, in Ref. it is assumed that 
2 < Xga < 4, which requires < 143 + 29x + 3x 2 - 1% < 120 + 18a;, which is not 
possible. Therefore the model is actually self-contradicting. However, Ref. fJJ has 
not stated this constraint on Xp- and it shall hence be ignored. This gives 4+3=7 
different models. The constraints on the X-charges of the right-handed neutrinos 
are < X^r < -X w < X L 2 < X w . 

[T2*] furthermore considers a model that was not dealt with in Ref. [11], aimed at 
the adiabatic Mikheev-Smirnow-Wolfenstein effect [HHl EH] for the solar neutrinos. 
The authors demand < X^2 < — X^s = X^i and the following constraints on the 
X-charges of the right-handed neutrinos: < X^r, X L 2 < — X^j < X L i < Xj^. 
Specifically they use X L i = |, X L 2 = |, Xp- = X-p- — 4. This results in z = 1 and 
X^r = Xp- ; furthermore the authors demand | < Xp- < y, which is fulfilled. So 
there are 4+3=7 models. 

[T5] indirectly states that a; = 3,y = 0,2; = and gives an explicit charged lep- 
ton mass matrix and a right-handed neutrino Majorana mass matrix. The latter 
depends on a parameter labeled X^, which is given as — |, in order to accomo- 
date the non-adiabatic MSW effect. This completely fixes the X L i, and one has 

y 13 y lj_ y 3 

N 1 ~ 2 '"^Af 2 ~~ 2 '^TV 3 ~~ 2- 

[Ti] states an explicit X-charge assignment, with right-handed neutrino charges 
X^pr = X^2 = X^3 = 0. Note that totally uncharged fields may give rise to 
dangerous tadpoles, e.g. see Refs. jHSl EH ESI EH| • 

[To] states an explicit X-charge assignment, with right-handed neutrino charges 
Xjpr = 8, Xjj2 = 4, X^3 = 0. These X-charges are the only difference to the model 
in Ref. ^1] . This model and the previous one are the only ones that have the quark 
mass matrices of Eq. (|3.51j) . 

[TH] explicitly demands conserved R p , however without checking whether this is a 
consequence of the given X-charge assignments. The authors also do not constrain 
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the X-charges by anomaly cancellation, which gives them more freedom. We recon- 
sider their model taking the Green-Schwarz mechanism into account. The authors 
indirectly state that z — l,y — and give several X-charges parameterized by 
the integers m,n,p,r. Then these parameters are fixed to be p = r = 5, n = 
(so-called "best-fit assignment"). One then distinguishes m = or 1 ("anarchical 
type" or "lopsided type", respectively), x and X L 3 are left undetermined. The 
authors thus have 2x4 = 8 models with X L 3 unconstrained, so one should al- 
ways be able to choose X L 3 such that R p is conserved. In the end of the pa- 
per the authors explicitly emphasize that tan/3 ~ 3 =3- x = 3, together with 
2Xq3 — Xjp — Xp- = 1. This seems to give a way to fix X L 3] however with Table 1 
one gets 2Xq3 — Xjj^ — Xp- = ^p, independent of X L 3 (note that one actually 
needs to have m — 1, otherwise the model is inconsistent). So we had to leave X^3 
unfixed. Although explicitly including right-handed neutrinos, no corresponding 
X-charges are stated. 



30 



by 


a; 


y 


z 




Xl2 




Compatible 
with Exp. 


O.K. if ... 


m 











X L 3 


X L 3 


1 -3n, 
ne 7L 


no if —3, . . . , 1 
no if n= —3, —2, —1 
no if n= 0, 1 


^ = 
A", ^ = 




2 








X L 3+1 


X L s-l 


3 — 3n, 
ne 7L 


no* if -2, . . . , 2 
no if n= —2, —1 
no if n= 0, 1, 2 


^ = 
A", * = 




2 








X L 3+1 


X L3 +5 


1 -3n, 
ne 7L 


no if —2, . . . , 2 
no if n= —2, —1 
no if n= 0, 1, 2 


^ = 
A", ^ = 


m 


2 








-12 


-13 


55 


no 


A" = 




2 








23 


22 


-78 


yes 















32 
5 


22 
5 


7 

5 


(no) h . . 


A' = 







-6 





26 
5 


16 

5 


1 
5 


no 


A' = 



continued next page 



31 



by 


X 


y 


z 




Xl-2 




Compatible 
with Exp. 


O.K. if ... 







-6 





62 
T 


48 
~ IT 


33 


no 


A, A' = 




1 


-6 





1321 
105 


989 
~ 105 


674 
_ 105 


no 


A, A' = 




2 


-6 





38 
~3~ 


28 


19 


no 


A, A' = 













33 
~ 


g 
5 


23 
IT 


(no) h . . 


A' = 







-6 





27 
~ 


2 

5 


17 
T 


no 


A' = 







-6 





5 


5 


5 


no 


A' = 













9Q 

5 


1 Q 

5 


4 
5 


(nojh.o. 


A' = 







-6 





zo 
5 


5 


2 

_ 5 


no 


A' = 







-6 





5 


5 


OtJ 

5 


no 


A, A' = 




1 


-6 





1258 
105 


1052 
105 


737 
_ 105 


no 


A, A' = 




2 


-6 





181 
"15" 


149 

~~ "15" 


104 
_ ~L5~ 


no 


A, A' = 













6 


1 


4 


(nojh.o. 


A, A' = 







-6 





24 
~5~ 


1 

~~ 5 


14 

T 


no 


A = 



continued next page 



32 



by 


X 


y 


z 




X.JJ2 




Compatible 
with Exp. 


O.K. if ... 







-6 





58 

IT 


37 


52 

~ IT 


no 


A, A' = 


■ j ■ 

i 


3 





1 


-8 


-8 


-4 


(no) h . . 


M/ = 










1 


-4 


-7 


9 

~2 


(no) h . . 


M/ = 




3 





1 


j 


-4 


_ 2 


no 


A = 


1 u ■ 











-8 


-8 


-8 


no 


^ = 


■ ^ ■ 








1 


X L z +2 




OtJO ^ ^ ^t/o 

105 — ••• — 105 


yes 




13 











-6 


-3 


-6 


no 


A', ^ = 













-5 


-5 


-5 


no 


M/ = 













-4 


-7 


-4 


no 


^ = 













-3 


-9 


-3 


no 


A, A', \I> = 













-10 


-4 


-1 


no 


A', \& = 













-10 


-4 


-10 


no 


^ = 













-9 


-6 


-9 


no 


\I> = 



continued next page 



33 



by 


X 


y 


z 


X L i 


Xl-2 




Compatible 
with Exp. 


O.K. if ... 













-8 


-8 


-8 


no 


* = 













-7 


-10 


-7 


no 


* = 




2 








-7 


-3 


-8 


no 


^ = 




2 








-6 


-5 


-7 


no* 


^ = 




2 








-5 


-7 


-6 


no* 


^ = 




2 








-4 


-9 


-5 


no 


^ = 




2 








-9 


-8 


-10 


no* 


^ = 




2 








-8 


-10 


-9 


no* 


^ = 










1 


-6 


-3 


-2 


no 


A, A', ^ = 










1 


-4 


-6 


-10 


no 


^ = 




1 





1 


-3 


-5 


-9 


no 


^ = 










1 


-10 


-1 


-9 


no 


A, ^ = 










1 


-8 


-6 


-6 


no 


^ = 



continued next page 



34 



by 


X 


y 


z 


X L i 


Xl2 




Compatible 
with Exp. 


O.K. if ... 




1 





1 


-8 


-4 


-5 


no 


A, ^ = 




1 





1 


-6 


-9 


-2 


no 


M/ = 




2 





1 


-8 


-2 


-4 


no 


A, ^ = 




1 





1 


-10 


-7 


-9 


no* 


M/ = 




2 





1 


-10 


-5 


-8 


no* 


^ = 




2 





1 


-8 


-10 


-5 


no* 


^ = 


El 


1 





1 


1 -LD 

30 


113 
no 

30 


113 

30 


yes 




i 











1 Q 

5 


1 Q 

5 


1 Q 

5 


no 


A, A', \1> = 


[HI1I2I 











X L 3 + 3 


X L 3 + 1 


/ 41 


1 

yes / no 


A, A' = 




1 








X L z + 3 


X L Z + 1 


/ 239 
•••/ ~~ 105 


1 

yes / no 


A, A' = 




2 








X L z + 3 


X L 3 + 1 


/ 11 


1 

yes / no 


A, A' = 




3 








X L z + 3 


X L 3 + 1 


/ 7 

•■■/-I 


1 

yes J no 


A, A' = 










1 


X L z + 3 


X L 3 + 1 


/ 73 
•••/ _ 35 


1 

yes / no 





continued next page 



35 



by 


X 


y 


z 


X L i 


^L 2 




Compatible 
with Exp. 


O.K. if ... 




1 





1 


X L3 + 3 


X L 3 + 1 


/ 33 
•••/ 20 


/ 

yes 1 no 






2 





1 


X L s + 3 


X L 3 + 1 


/ 11 


1 

yes 1 no 






3 





1 


X L3 + 3 


X L 3 + 1 


/ 4 


1 

yes / no 









-6 





X L z + 3 


X L 3 + 1 


/ 59 
-/ - 15 


I 

yes / no 


A, A' = 




1 


-6 





X L 3 + 3 


X L 3 + 1 


/ 347 
•••/ _ T05 


1 

yes / no 


A, A' = 




2 


-6 





X L 3 + 3 


X L 3 + 1 


/ 41 

■■■/-* 


1 

yes / no 


A, A' = 







-6 


1 


X L 3 + 3 


X L 3 + 1 


/ i no 
"7 35 


1 

yes / no 






1 


-6 


1 


X L s + 3 


X L 3 + 1 


•••/ 20 


I 

yes / no 






2 


-6 


1 


X L 3 + 3 


X L 3 + 1 


"■'/ 45 


1 

yes / no 




EH EE] 








1 


5 


2 


-2 


yes 






1 





1 


5 


2 


-2 


yes 






2 





1 


5 


2 


-2 


yes 






3 





1 


5 


2 


-2 


yes 





continued next page 
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by 


X 




z 


X L i 






Compatible 
with Exp. 


O.K. if ... 







-6 


1 


5 


2 


-2 


yes 






1 


-6 


1 


5 


2 


-2 


yes 






2 


-6 


1 


5 


2 


-2 


yes 




H2! 








1 


9 
2 


3 
2 


9 

~2 


yes 






1 





1 


9 
2 


3 
2 


9 

_ 2 


yes 






2 





1 


Q 

y 
2 


Q 
O 

2 


Q 

y 
2 


yes 






3 





1 


Q 

y 
2 


Q 
O 

2 


Q 

y 
2 


yes 









-6 


1 


Q 

y 
2 


Q 
O 

2 


Q 

y 
2 


yes 






1 


-6 


1 


Q 

y 
2 


Q 
O 

2 


Q 

y 
2 


yes 






2 


-6 


1 


9 
2 


3 
2 


g 

_ 2 


yes 




E3! 


3 








79 
30 


1\ 
~ 30 


_ 30 


(no) h . . 


* = 


■ j j ■ 

M 





1 














no 









1 














no 





continued next page 
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by 


X 


y 


z 


X L i 


Xl2 


Xls 


Compatible 
with Exp. 


O.K. if ... 










1 


X L 3 


X L 3 




yes 






1 





1 


X L 3 


X L 3 




yes 






2 





1 


X L 3 


X L 3 




yes 






3 





1 


X L 3 


X L 3 




yes 












1 


X L 3 + 1 


X L 3 




yes 






1 





1 


X L 3 + 1 


X L 3 




yes 






2 





1 


X L 3 + 1 


X L 3 




yes 






3 





1 


X L 3 + 1 


X L 3 




yes 





Table A: In the table above we present the catalogue of Froggatt-Nielsen models of the 
type defined in Section 1, parametrized by x, y, z, X L i, X L 2, X L 3. For the entries "yes", 
"no", etc. and \& = ; A" = etc. see Section 6. Models where the "no" is labeled by a 
* are more in accord with the bounds if rh > 250 GeV, one gets a (no)h.o. instead. 
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by 


X H T> , X H U 


X 10 i 


^5* 


v 

AT 4 


e 


Compatible with 
Constraints? 


[16J 


O . o 

2g 3 , -2g 3 


O i O i 

3 + q 3 ,2 + q 3 ,q 3 


1 i ,J J J 

1 + d 3 , cf 3 , d 3 




0.22 


no 


1 *7 I 

mi 


n n 
(J, (J 


o 1 n 
2, 1, U 


inn 
1, U ; (J 


0, 1,2 


1 


no 




n n 
U, U 


o 1 n 
Z, 1, U 


inn 
1, U, U 


n n o 

0,0,2 


l 

V300 


no 




n n 

u, u 


o 1 n 

z, 1, U 


inn 

1, 0, 


1 1 o 

1, 1,2 


1 

Vmj 


no 




n n 
(J, (J 


z, 1, U 


O 1 1 

Z, 1, 1 


n n i 

0,0, 1 


i 

Vmi 


no 




n n 
U, U 


o 1 n 
Z, 1, U 


inn 
1, U, U 


n n 1 

0,0, 1 


i 

V300 


no 


I 1 o I 

[18J 


u, u 


o o n 

3, 2, 


o n n 

2, 0, 




n oo 

U.zz 


no 




n n 
(J, (J 


o o r» 

6 } z, (J 


O 1 1 

0, 1, 1 




n oo 

U.zz 


no 




n n 
U, U 


o o n 
o, Z, U 


A O O 

4, z, z 




n oo 
U.ZZ 


no 


1 1 n I 

[19J 


n n 

u, u 


o o n 

3, 2, 


o n n 

2, 0, 


o n n 

2,0,0 


n oo 

U.zz 


no 


I on I 

[20J 


n n 
(J, (J 


2, 1, 


1 n n 

1, U ; (J 


o 1 n 

2, 1,0 


n n*7 

U.U7 


no 




n n 
U, U 


2, 1, U 


O 1 1 

Z, 1, 1 


o i n 

2, 1,0 


n ri7 
U.U f 


no 




n n 

u, u 


2, 1, 


inn 

1, 0, 


n n n 

0, 0, 


n n^7 

0.07 


no 




n n 

u, u 


2, 1, 


2, 1, 1 


n n n 
(J, (J, (J 


n n*7 

0.07 


no 




U, U 


3, 1, U 


O 1 1 

2, -L, 1 


o i n 

2, 1,0 


U.U f 


no 




U, U 


o, 1, U 


O 1 1 

z, 1, 1 


n n n 
U, U, U 


U.U f 


no 


m 


0,0 


2,1,0 


1,1,1 


0,0,0 


U.U5 


no 



Table B: The models in this table do not lead to a fermionic mass spectrum as in 
Eqs. \1.V\1.1\) . Note that all models are SU(5) invariant. 



B Examining the Mass Matrices 

To demonstrate the validity of Eqs. ()3.49H3.52l . note that the CKM matrix is nearly a 
unit-matrix, so that one should have U^ Ul ^ ~ U^ Dl ^\ i.e. these two matrices should 
have "the same Cabibbo structure" (for a more detailed analysis see Ref. JH]). It is 
thus reasonable that G^G^ « G^G^\ This is indeed the case for the matrices 
presented, e.g. Eq. (j3.51j) gives 
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Furthermore the £r'"") have the correct eigenvalues 27 and hence reproduce the correct 
masses. For example, in Eq. (|3.49jl has the characteristic polynomial 

-det^e 12 + g^% 3 - g^ 23 g^\ 2 ) e 4 

(9 {U) n 9 (U \ 3 ~ 9 (U \ 3 9 {U \.) e* + (<?% g^ 22 - «,% 9 



21; 



A 



Assuming that 



-(y% 6 8 + 0% e 4 + g^ m ) A 2 + A 3 = 0. 

(B.3) 



det g (u) « 1 , 

9 22 y 33 y 23 y 32 ~ 1 



a {u) a (v) _ a {u) a {u) ~ I 
y 33 9 11 y 31 9 13 ~ 1 > 

y% 9 iU \ 2 ~ 9 iU) n 9 iU) 21 - 1 , 

9 { %s - 9 { % 2 - 2 (C/ 11 - 1 , (B.4) 



this reduces to 



A 3 -(e 8 + e 4 + l) A 2 + (e 4 + e 8 + e 12 ) A - e 12 ^ , (B.5) 
which can be factored as 

(A-l)(A-e 4 )(A-6 8 )^0, (B.6) 

thus reproducing the correct ratio of the eigenvalues, see Eq. (1.5). Alternatively, in the 
expression for the characteristic polynomial one can also neglect higher order terms and 



27 As the matrix G^"") ■ G^ - ") is Hermitian, it can be diagonalized with its eigenvalues on the diagonal. 
Thus the matrix D^'"' which one gets when bi-unitarily diagonalizing G^'") contains the square roots 

of the eigenvalues of G^'") ■ (times a phase) on its diagonal. This should not be confused with 

the strictly speaking wrong statement that £)(•••) contains the eigenvalues of G^'"\ As a demonstration, 
consider the following matrix, its coefficients having been randomly generated with Mathematical® , their 



absolute values being in the interval 



0.8614-2.9121 i -1.2754 + 0.1818 i 0.6656 - 2.2503 i 
.V I [ -1.2074- 1.3980 i 1.2850 - 0.3739 i 0.5473 + 0.8524 % 
-0.2591 + 0.5662 i -0.4475 + 0.0637 i -0.6963 + 1.5219 i 

The absolute values of the eigenvalues of A4 are 2.6456, 2.5975, 1.7199. However, the absolute values 
of the square roots of the eigenvalues of are 4.3262, 2.4631, 1.1091. On the other hand, if the 

entries of Ai are e-suppressed, e.g. with the exponents 




one gets 2.3586, 0.0813, 0.0051 and 2.3701, 0.0819, 0.0050, respectively, almost the same. 
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only demand 



o iu) « 1 
y 33 ~ x ' 

o (*0 „{u) _ n {u) (to ~ i 
y 22 y 33 y 23 y 32 ~ x ) 

detc/ (c/) ^l, (B.7) 

so that 

A 3 - A 2 + e 4 A - e 12 = . (B.8) 

Using e = 0.22 we thus obtain Ai = e 8 - 00212 ) \ 2 = e 3 -" ) A 3 = e - 001 , which is also in good 
agreement. 

To numerically demonstrate that G (C/) ,G (Z?) of Eq. give the correct CKM- 

matrix, random Oil) coefficients were generated with Mathematica® . With Eq. ()3.49j) 
we obtained 

= 
= 

Using e = 0.22, this leads to 

m b :m s :m d = -1.98 : -0.68 e 2 : 1.13 e 4 , (B.10) 



1.41 e 8 -0.95 e 5 -1.73 e 3 \ 
0.94 e 7 -1.19 e 4 1.12 e 2 , 
1.35 e 5 2.10 1.52 e 2 / 

1.22 e 4 -1.45 e 3 -1.76 e 3 \ 

0.69 e 3 -0.94 e 2 0.60 e 2 . (B.9) 
-1.02 e 0.68 -1.97 / 



m t :m c :m u = -1.51 : 0.33 e 4 : 2.56 e 8 , (B.ll) 

and 

/ -0.98 -0.62 e 3.66 e 3 \ 

jjCKM = Q 62 £ _ Q Q8 _ Q 8Q e 2 ^ (B.12) 

\ 4.18 e 3 -0.77 e 2 0.99 / 

which is the correct order of magnitude. 

Finally for completeness it should be mentioned that Eqs. ()3.49j) and (|3.51j) can also be 
generated from models with a vector-like pair of Froggatt-Nielsen fields A, B, see Ref. jTj. 
A vector-like pair of Froggatt-Nielsen fields can furthermore generate phenomenologically 
viable mass matrices which cannot be generated from models with only one Froggatt- 
Nielsen field: e.g. Xqi = Xjjt = —8, Xq2 = Xjp = 2, Xq3 = Xjjj = X H u = and 
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[A) = (B), see Ref. |2Ej, give : 



£ 16 £ 6 £ S 



oc e 6 e 4 e 2 . (B.13) 



e 8 e 2 e° 



This leads to a characteristic polynomial of approximately the form 

A 3 - (l + e 4 + e 16N )A 2 + (V + e 12 + e 16 + e 20 ) A - (> + e 16 + e 20 ) = 0, (B.I 1) 



with e = 0.22 giving 

A 1 = l, A 2 = e 4 - 00156 , A 3 = e 7 - 99689 (B.15) 

and thus 

m u : m c : m t = e 8 : e 4 : 1, (B.16) 

as required. 
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